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We numerically investigate localization properties of electronic states in a static model of 
poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers with realistic parameters obtained by 
quantum-chemical calculation. The randomness in the on-site energies caused by the electron- 
phonon coupling are completely correlated to the off-diagonal parts. In the single electron model, 
the effect of the hydrogen-bond stretchings, the twist angles between the base pairs and the finite 
system size effects on the energy dependence of the localization length and on the Lyapunov expo- 
nent are given. The localization length is reduced by the influence of the fluctuations in the hydrogen 
bond stretchings. It is also shown that the helical twist angle affects the localization length in the 
poly(dG)-poly(dC) DNA polymer more strongly than in the poly(dA)-poly(dT) one. Furthermore, 
we show resonance structures in the energy dependence of the localization length when the system 
size is relatively small. 
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> 

OO 

o 

o 

-i— > 



T3 

c 

o 
o 



X 



I. INTRODUCTION 



The recent development of the nanoscale fabrication let 
us expect the utilization of the DNA wire as a molecular 
device 0, and the realization of DNA computing Q . 
DNA is believed to form an effectively one-dimensional 
molecular wire, which is highly promising for diverse ap- 
plications. Actually, the modern development of physico- 
chemical experimental techniques enables us to measure 
directly DNA electrical transport phenomena even in sin- 
gle molecules 0, 0- Moreover, several groups have re- 
cently performed numerical investigations of localization 
properties of DNA electronic states based on realistic 
DNA sequences HH0 ■ 

However, DNA transport properties still remain a con- 
troversial topic, mainly due to tremendous difficulties 
in setting up the proper experimental environment and 
the complexity of the molecule itself. Specifically, a dis- 
tinctive feature of biological polymers is the complicated 
composition of their elementary subunits, and an appar- 
ent ability of their structures to support long-living non- 
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linear excitations. Although a number of theoretical ex- 
planations for DNA charge transfer/transport phenom- 
ena have been suggested on the basis of standard solid- 
state- physical approaches, like polarons, solitons, elec- 
trons or holes [Ilill[l3,[ll[ll[ll[ll[l3,the situation 
is still far from working out a unique, non-contradictory 
theoretical scheme. 

In their polaron-like model, Hennig and coworkers 
studied electron breather propagation along DNA ho- 
mopolynucleotide duplexes, i.e. in both poly(dG)- 
poly(dC) and poly(dA)-poly(dT) DNA polymers ||, and, 
for this purpose, estimated electron-vibration coupling 
strength in DNA using semiempirical quantum-chemistry 
[Til Il7j . Chang et al. have also considered a possible 
mechanism to explain the phenomena of DNA charge 
transfer 0. The charge coupling with DNA structural 
deformations can create a polaron and thus promote a 
localized state H[I3. As a result, the moving electron 
breather may contribute to highly efficient long-range 
conductivity. Recent experiments seem to support the 
polaron mechanism for the electronic transport in DNA 
polymers [Isj . 

In the present paper, we investigate localization prop- 
erties of electronic states in a stochastic bond vibra- 
tion model of poly(dG)-poly(dC) and poly(dA)-poly(dT) 
DNA polymers by adopting the model by Hennig et al. 
PI. Here we assume that the disorder is caused by DNA 
vibrational modes, and via electron-vibrational coupling 
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it influences the charge transfer/transport along DNA 
duplexes. Moreover, we discuss the difference in charge 
localization properties between poly(dG)-poly(dC) and 
poly(dA)-poly(dT) DNA polymers, as well as the pecu- 
liarities of the mixed (AT/GC) model. 

The outline of the present paper is as follows. In the 
next section we introduce the DNA model to investigate 
subsequently the electronic states and their coupling to 
DNA vibrational modes. In the Sec. 3 we present nu- 
merical results concerning the influence of changes of the 
hydrogen-bond stretchings and twist angles as well as the 
effects of finite system size on the localization properties. 
The last section contains a summary and discussion. 

II. MODEL 

This section reformulates the model Hamiltonian for 
DNA duplexes, as introduced by Hennig et al. Q 
into a one-electron tight-binding adiabatic Hamiltonian 
with structural disorder, taking into account a nearest- 
neighbor electron hopping term. The latter model is ac- 
tually akin to that proposed earlier by Dunlap-Kundu- 
Philips (DKP). Specifically, the general DKP model con- 
tains disorder induced by the proper vibrational modes, 
which shows up in both the diagonal and the off-diagonal 
parts of the electron Hamiltonian [lj| . Generally, in the 
DKP model the transport can arise from a set of measure- 
zero unscattered states at a particular energy derived 
from the resonant transmission formula |l9l |20| , and the 
number of the unscattered states can be estimated by the 
resonance width in the system size. 

The Hamiltonian for the electronic part in our DNA 
model is given by 

H e i = ^ E n C\C n 

n 

- £vwi(ctc n+1 + c n+1 ct), (1) 

n 

where C n and are creation and annihilation operators 
of an electron at the site n. The on-site energies E n are 
represented as 

E n = E + kr n , (2) 

where Eq is a constant and r„ denotes the fluctuation 
caused by the coupling with the transversal Watson-Crick 
H-bonding stretching vibration. 

The transfer integral V nn +\ depends on the three- 
dimensional stacking-distance d nn +i between adjacent 
bases labeled by n and n + 1, along each strand and 
is given as, 

V nn+ i = V (l - ad nn+ i) . (3) 

The parameters k and a, which describe the strength of 
the interaction between the electronic and spatial vari- 
able, have been previously calculated through quantum- 
chemical methods. Radial displacements bring about 



also a variation of the distances between neighboring 
bases along each strand d nn+ i. The first-order Taylor 
expansion around the equilibrium positions is given by 

R 

d n n+i = -^-(1 - cos8 )(r n + r n+1 ). (4) 
«o 

Rq represents the equilibrium radius of the helix, 9q is 
the equilibrium twist angle between base pairs, and £q 
the equilibrium distance between bases along a strand 
given by 

£o = (a 2 +4i?2 s in 2 (V2)) 1/2 , (5) 

with a being the distance between neighboring base pairs 
in the direction of the helix axis. A sketch of the geo- 
metrical parameters Ro,£q,9o, r n +\ and d nn+ i is given in 
Fig.l. 




FIG. 1: Sketch of the structure of the DNA model. The bases 
are represented by bullets and the geometrical parameters 
Ro, £o, do, r„+i and d nn+1 are indicated. 

We adopt realistic values of the parameters obtained 
from the semi-empirical quantum-chemical calculations; 
k = O^SWleVk- 1 }, a = 0.053835L4.- 1 ] for the cou- 
pling parameters of the poly(dA)-poly(dT) DNA poly- 
mer, k = -0.090325[eVA~ 1 ], a = 0.383333IA- 1 ] for 
the ones of the poly(dA)-poly(dT) DNA polymer. As 
to the other typical parameters for DNA molecules, we 
use: E = 0.1[eV],V = 0.1[eV], a = 3.4[A],i? = 10 [A] 
and O = 36°. 

Further, we consider {r n } as independent random vari- 
ables generated by uniform distribution with width (r n 6 
[— W, W]). Accordingly, fluctuations in both the on-site 
energies and the off-diagonal parts in the Hamiltonian 
(1) are mutually correlated because they are generated 
by the same random sequence r n . (See Fig. 2(c).) The 
typical value for W is W = 0.1 [A], which approximately 
corresponds to the variance in the hydrogen bond lengths 
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in Watson- Crick base pairs, as seen in X-ray diffraction 
experiments [21| . 

Figures 2(a) and (b) show a typical fluctuation pattern 
of E n and V nn +i for W = 0.1 in poly(dG)-poly(dC) and 
poly(dA)-poly(dT) DNA polymers. Although the fluc- 
tuation of the on-site energy is almost of the same or- 
der in both poly(dG)-poly(dC) and poly(dA)-poly(dT) 
DNA polymers, the fluctuation of the transfer integral 
Vnn+i m the poly(dG)-poly(dC) DNA polymers is much 
larger than that of the poly(dA)-poly(dT) DNA poly- 
mers. This property ought to reflect the difference in the 
electron localization nature between poly(dG)-poly(dC) 
and poly(dA)-poly(dT) DNA polymers. We can see the 
correlation between the sequences, E n , K m +i, in the 
parametric plots in Fig. 2(c). 

In addition to the DNA homopolymer duplexes, we in- 
vestigate the localization nature of the mixed sequence 
consisting of two types of the Watson-Crick pairs. Then, 
as a zero-order approximation, the electron-phonon cou- 
pling parameters for the mixed GC/AT stacks are taken 
here to be equal to the values obtained for poly(dG)- 
poly(dC) and poly(dA)-poly(dT) DNA polymers. 



III. NUMERICAL RESULTS 

The Schrodinger equation H e i\& >= E\Q > is written 
in the transfer matrix form, 

(r)=(¥ "¥)(:--,)• <« 

where a n is the amplitude of the electronic wavefunction 
$) = J2 n an \ n ) at the base pair site n. We use the 
localization length £ and/or Lyapunov exponent 7 cal- 
culated by the mapping (6) in order to characterize the 
exponential localization of the wave function. Originally 
the Lyapunov exponent (inverse localization length) is 
defined in the thermodynamic limit (N — > 00), however, 
in the present paper we use the following definition for 
the Lyapunov exponents of the electronic wave function 
for a large system size N [22, |2ij ■ 

l{ E,N)=C\E,N) ^ a ^^ a ^ 2 \ (7) 

We use appropriate initial conditions <xq = a\ = 1, and 
for large N(>> £) the localization length and Lyapunov 
exponent are independent of the boundary condition. 
The system size dependence of the localization length 
is investigated in Subsect.3.3 in the relation with the res- 
onance energy. The energy-dependent transmission coef- 
ficient T(E, N) of the system between metallic electrodes 
is given as T{E, N) = exp(— 2^N) and is related to Lan- 
dauer resistance via g = (1— TVT in unit of the quantum 
resistance h/2e 2 (~ 13[kflm\) |24| . 
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FIG. 2: (a)The on-site energy _E„[eV] and (b)transfer integral 
Vnn+i [e V] as a function of the base pair site n. The paramet- 
ric plot E n versus Vnn+i is shown in (c). W = 0.1 and the 
other parameters are given in text. The unit of the energy 
and the spatial length are [eV] and the number of nucleotide 
base pair [bp], respectively, throughout the present paper. 



A. Fluctuation effect of stretching 

Figure 3(a) shows the dependence of the localization 
length £ on the energy for the cases of W = 0.1, 0.2, 0.4 in 
the poly(dA)-poly(dT) DNA polymer. For the computa- 
tion of the energy dependence of the localization length 
we used a energy subdivision of 800 points in the en- 
ergy range. The extent of the localization depends on 
the Fermi energy of the electroads when we measure the 
transmission coefficient or conductance. It is apparent 
that the localization length is reduced as a result of the 
increase of the fluctuation strength W of the hydrogen 
bond stretching. The small peaks for the resonance en- 
ergy are suppressed in W > 0.2 due to the large fluctua- 
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tion in the r n . 

Figure 3(b) and (c) show the localization length as a 
function of the energy in the case of poly(dG)-poly(dC) 
and the mixed DNA polymers, respectively. In Fig. 3(b) 
the global behavior is almost the same as that illustrated 
in Fig. 3 (a) except for the shape of the energy depen- 
dence. The localization length in the poly(dA)-poly(dT) 
DNA polymer is globally larger than the one in poly(dG)- 
poly(dC) DNA polymer. The reason is simply because 
the fluctuation strength of V nn +i m poly(dG)-poly(dC) 
DNA polymers is larger than the one in poly(dA)- 
poly(dT) DNA polymers as seen in Fig. 2. 

In the mixed case, the localization length in the mixed 
sequence shows the intermediate behavior between the 
poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA poly- 
mers. The resonance peak structure in the energy depen- 
dence of the localization length is not as pronounced as 
in the case of poly(dA)-poly(dT) and poly(dG)-poly(dC) 
DNA polymers. 

As a result, more disorder in the H-bond stretching 
(higher W values) renders the peak structure in the elec- 
tronic states to disappear and makes the latter more lo- 
calized. The finite system effect will be discussed in Sub- 
sect. 3. 3. 



B. Twist angle effect 

In this subsection, we investigate how the localization 
length depends on the helical twist angle Oq of the DNA 
double helix. It is known that the allowed values of the 
angle range from ca. 27° till 45°. 

Figure 4 shows the energy dependence of the localiza- 
tion length for some twist angles in the DNA models. 
We conclude, that in poly(dA)-poly(dT) polymers the 
change of the angle does virtually not affect the local- 
ization length, however, in the poly(dG)-poly(dC) poly- 
mer it does have an impact. The difference is caused by 
the 9q dependence of the fluctuation of the transfer en- 
ergy V nn+ i which is smaller in the case of the poly(dG)- 
poly(dC) polymer than in its poly(dA)-poly(dT) coun- 
terpart as seen in Fig. 2(b). 

It follows that when the twist angle increases in the 
poly(dG)-poly(dC) DNA polymer the localization length 
becomes smaller (larger) in the low (high) energy region. 



C. System size effect on the resonance 

Figure 5 shows the system size dependence of the lo- 
calization length in the poly(dG)-poly(dC) polymer and 
the poly(dA)-poly(dT) polymer. Although the states for 
most of energies are exponentially localized, and have 
well-pronounced resonant maxima at discrete points cor- 
respond to the eigenenergies of each system when the 
system size is relatively small (N = 2 16 ). Electronic 
states whose energy is close to resonance peaks have a 
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FIG. 3: The Localization length as a function of energy for 
several values of W{= 0.1,0.15,0.2) in the (a) poly(dA)- 
poly(dT) DNA polymers, (b) poly(dG)-poly(dC) DNA poly- 
mers and (c) mixed DNA polymers. The system size is 
TV = 2 22 and 8q — 36°. 



tendency to be extended in the sense that their local- 
ization length is larger than the system size. For the 
larger system size the energy dependence of the localiza- 
tion length converges. 

In our definition of the localization length, when the 
system size N is smaller than the localization length, the 
sharp energy dependence of the resonance peaks appears, 
and the structure of the resonance depends on the ran- 
domness in r n of each sample. Note that the resonance 
structure can be washed out if we take the ensemble av- 
erage < £ > with regard to different configurations of r n 
even for W — 0.1. 

As a result we can say the transfer and transport of 
the electron sharply depend on the Fermi energy of the 
electroades as the system size is relatively small. The 
extended states caused by the resonance can contribute 
to the transfer of the electron. 
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FIG. 4: The localization length as a function of energy for 
several values of the twist angle 9 {= 0° , 27°, 36° , 45°, 60°) 
in the (a) poly(dA)-poly(dT) DNA polymers, (b) poly(dG)- 
poly(dC) DNA polymers and (c) mixed DNA polymers. The 
system size is N = 2 22 and W = 0.1. 



IV. SUMMARY AND DISCUSSION 

We have numerically investigated localization proper- 
ties of electronic states in an adiabatic polaron model of 
poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA poly- 
mers with realistic parameters obtained using semi- 
empirical quantum-chemical calculations. 

Now, we are in a position to compare the localization 
properties of the poly(dG)-poly(dC), poly(dA)-poly(dT) 
DNA polymers and the mixed case. Figure 6(a) and (b) 
show the localization length and Lyapunov exponent (in- 
verse of the localization length) in the three types of the 
polymers with W = 0.1. In the low energy region, the 
localization length in the poly(dG)-poly(dC) DNA poly- 
mer is larger than that in the poly(dA)-poly(dT) DNA 
polymer. On the other hand, it is known from the experi- 
ments that poly(dG)-poly(dC) oligomer is a better semi- 
conductor than poly(dA)-poly(dT) oligomer [lSj. The 
difference is caused by the dynamical effect which is ne- 
glected in this paper and the system size effect. Indeed, 
the localization length of the DNA polymers is larger 
than £ > 2000 [bp] in almost all the energy band for all 
models, that is much larger than the system size of the 
oligomer used in the experiments. As seen in Fig. 5 in the 
smaller system size the resonance peaks become complex. 
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FIG. 6: Comparison (a)localization length and (b)Lyapunov 
exponent in poly(dG)-poly(dC), poly(dA)-poly(dT) and the 
mixed DNA polymers. W = 0.1, 9 = 36° and N = 2 22 . 
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